Hofava-Lifshitz cosmology with generalized Chaplygin gas 
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We investigate cosmological scenarios of generalized Chaplygin gas in a universe governed by 
Hofava-Lifshitz gravity. We consider both the detailed and non-detailed balance versions of the 
gravitational background, and we include the baryonic matter and radiation sectors. We use obser- 
vational data from Type la Supernovae (SNIa), Baryon Acoustic Oscillations (BAO), and Cosmic 
Microwave Background (CMB), along with requirements of Big Bang Nucleosynthesis (BBN), to 
constrain the parameters of the model, and we provide the corresponding likelihood contours. We 
deduce that the present scenario is compatible with observations. Additionally, examining the evolu- 
tion of the total equation-of-state parameter, we find in a unified way the succession of the radiation, 
matter, and dark energy epochs, consistently with the thermal history of the universe. 
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I. INTRODUCTION 

Recently Hof ava proposed a power-counting renormal- 
izable theory with consistent ultra-violet (UV) behav- 
ior Although presenting an infrared (IR) fixed 
point, namely General Relativity, in the UV the the- 
ory exhibits an anisotropic, Lifshitz scaling between time 
and space. Due to these novel features, there has been 
a large amount of effort in examining and extending 
the properties of the theory itself |5l— [30| - Additionally, 
application of Hofava-Lifshitz gravity as a cosmological 
framework gives rise to Hofava-Lifshitz cosmology, which 
proves to lead to interesting behavior [3l|, HH . In partic- 
ular, one can examine specific solution subclasses [33l — 
|40| . the phase-space behavior (4lH45j . the gravitational 
wave production [46H49|. the p erturbation spectrum (50r - 
l5al . the matter bounce |59H63l ] , the black hole properties 
[64H72| , the dark energy phenomenology |73H76| , the ob- 
servational constraints on the parameters of the theory 
[77l l78j , the astrophysical phenomenology (79H8r| , the 
thermodynamic properties |87H9l| etc. However, despite 
this extended research, there are still many ambiguities 
if Hofava-Lifshitz gravity is reliable and capable of a suc- 
cessful description of the gravitational background of our 
world, as well as of the cosmological behavior of the uni- 
verse @, H, EH, [H, [93-IH . 

Although the foundations and the possible conceptual 
and phcnomcnological problems of Hofava-Lifshitz grav- 
ity and its associated cosmology are still an open issue, 
it is worth investigating different cosmological scenarios 
in this gravitational background. In this regard, it would 
be interesting to examine cosmological scenarios where 



apart from the gravitati onal sector there exists a gener- 
alized Chaplygin gas [9ril - ll03| . 

The generalized Chaplygin gas (GCG) is an alterna- 
tive to the Quintessence model which had attracted great 
interest in recent times. The scenario can explain the ac- 
celeration of the universe via an exotic equation of state, 
which mimics a pressureless fluid at the early stages of 
evolution of the Universe, and a cosmological constant 
at late times. It is therefore interesting to consider the 
GCG scenario as a unified description for dark matter 
and dark energy (98|. Th e ba c kgro und evolution fits 
well the observational data [lOOl . Il0l| , however the cos- 
mological behavior is indistinguishable from that of the 
ACDM scenario while fitting with the structure forma- 
tion data as well as w ith data f rom Cosmic Microwave 
Background radiation |l04l . 1 1051 ] . Additionally, the sce- 
nario is plagued by the presence of instabilities as well as 
oscillations which are not observed in the matter power 
spectrum. Although this is a serious drawback of GCG 
scenario it is not the final verdict for its fate as a model for 
the unified description of dark matte r and dark energy. 
As it was shown by Reis et al. in |106| , allowing for small 
entropy perturbations can eliminate instabilities and os- 
cillations in the matter power spectrum, even in the lin- 
ear regime, for a region of the parameter space where the 
GCG model behaves quite differently from the ACDM 
one. Furthermore, as it was shown by Avelino et al. in 
1071 ] , in the GCG scenario the transition from dark mat- 
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tcr to dark energy behavior is never s mooth, an d hence 
the linear theory which was used by pM [lOi] to rule 
out GCG as a unified model, may break down late in the 
matter-dominated era, even on large cosmological scales. 
Therefore, nonlinear effects should be necessarily taken 
into account when confronting cosmological observations. 
Moreover, the addition of baryons in the GCG scenario 
can al so im prove the behavior of the matter power spec- 
trum 108]. In summary, the GCG as a unified model 
of dark matter and dark energy is not completely ruled 
out and it deserves further investigations. Finally, we 
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mention that GCG in the presence of cold dark matter 
as well as baryonic matter (where GCG acts as a normal 
dark energy candidate) is one of the most well-fit sce- 
narios with cosmological observations, amongs t all the 
exotic models that have been considered so far |l Q9j ] . 

In the present work we construct the cosmology of 
a generalized Chaplygin Gas in a universe governed 
by Hofava-Lifshitz gravity. Additionally, we use obser- 
vational data from Type la Supernoyae (SNIa) |l!0j |. 
Baryon Acoustic Oscillations ( BAO ) and Cosmic 

Microwave Background (CMB) |l!2j |. together with the 
Big Bang Nucleosynthesis (BBN) conditions, to constrain 
the various parameters of the model. Furthermore, in 
order to be general we perform the analysis both with 
and without the detailed-balance condition of the gravi- 
tational sector. 

The manuscript is organized as follows: In sectionllTlwc 
present Hofava-Lifshitz cosmology in both its detailed- 
balance and beyond-detailed-balance version. In section 
IIIII we construct the scenario of a generalized Chaply- 
gin Gas in Hofava-Lifshitz gravitational background and 
we extract the cosmological equations. In section HVl we 
use observational data in order to constrain the various 
parameters of the scenario. In section fVl we discuss the 
cosmological implications, focusing on the evolution of 
the total equation-of-state parameter of the universe and 
of the expansion rate. Finally, section IVII is devoted to 
the summary of our results. 



II. HORAVA-LIFSHITZ COSMOLOGY 

In this section we briefly review the scenario where 
the cosmolo gica l evolution is governed by Hofava-Lifshitz 
gravity [UHIi]- The dynamical variables are the lapse 
and shift functions, N and iVj respectively, and the spa- 
tial metric gij (roman letters indicate spatial indices). In 
terms of these fields the full metric is written as: 



ds 2 = -N 2 dt 2 + g ij (dx i + N l dt){dx 3 + N j dt), (1) 

where indices are raised and lowered using <?.y . The 

scaling transformation of the coordinates reads: t — > 
l 3 t and x % -> lx % . 



the full action of Hofava-Lifshitz gravity is given by 



S 



J dtd 3 x^N ^(KijK ij - XK 2 ) 



8(3A-1) 



2w 2 ^g 



1 - 4A 



-R 2 + AR- 3A' 



where 



Ki 



2N 



is the extrinsic curvature and 



e ijk 



(2) 



(3) 



(4) 



the Cotton tensor, and the covariant derivatives are de- 
fined with respect to the spatial metric <7y. e^ k is the 
totally antisymmetric unit tensor, A is a dimcnsionlcss 
constant and the variables k, w and fx are constants with 
mass dimensions —1, and 1, respectively. Finally, we 
mention that in action @ we have already performed 
the usual analytic continuation of the parameters p and 
w of the original version of Hofava-Lifshitz gravity, since 
such a procedure is required in order to obtain a realistic 
cosmology [HI, H3, IH, H3] (although it could fatally affect 
the gravitational theory itself). Therefore, in the present 
work A is a positive constant, which as usual is related 
to the cosmological constant in the IR limit. 

In order to add the matter component we follow the hy- 
drodynamical approach (which is most suitable for phe- 
nomenological analysis) of adding a cosmological stress- 
energy tensor to the gravitational field equations, by de- 
manding to recover the usual general relativity formula- 
tion in the low-energy limit [H, |4l|, [7f| . Thus, this matter- 
tensor is a hydrodynamical approximation with p m and 
p m (or p rn and w m ) as parameters. Note that this p m 
is the total matter energy density, that is it accounts for 
both the baryonic pb as well as the dark matter pdm- Sim- 
ilarly, one can additionally include the standard-modcl- 
radiation component (corresponding to photons and neu- 
trinos) , with the additional parameters p r and p r (or p r 
and w r ). 

In order to investigate cosmological frameworks, we 
impose the projectability condition Q and we use an 
FRW metric 



Detailed Balance 



N = l, ft,- = a 2 (f)7« , N* = 0, 



(5) 



The gravitational action is decomposed into a kinetic 
and a potential part as S g = J dtd 3 Xy/gN(£K+Cv)- The 
assumption of detailed balance 0] reduces the possible 
terms in the Lagrangjan, and it allows for a quantum 
inheritance principle since the (D + l)-dimensional 
theory acquires the renormalization properties of the D- 
dimensional one. Under the detailed balance condition 



with 



dx l dx 3 = 



dr 2 



1 - Kr 2 



r 2 dVi\ , 



(6) 



where K <,=, > corresponding to open, flat, and 
closed universe respectively (we have adopted the con- 
vention of taking the scale factor a(t) to be dimcnsionlcss 
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and the curvature constant K to have mass dimension 2). 
By varying TV and g^j, we extract the Friedmann equa- 
tions: 



H 2 = 



6(3A - 1) 

K 2 

6(3A-1) 

k a h 2 AK 
8(3A- l) 2 a 2 



3^V 2 X 2 3/^VA 2 
8(3A- l)a 4 + 8(3A - 1) 



(7) 



H 



-H 2 



4(3A- 



K 2 fi 2 K 2 



3^VA 2 



4(3A - 1 

n A p 2 AK 
16(3A- l) 2 a 2 



8(3A-l)a 4 8(3A-1) 



(8) 



where H = ^ is the Hubble parameter. 

The term proportional to a~ 4 is the usual "dark radia- 
tion term" , present in Hofava-Lifshitz cosmology [3ll. |32| . 
while the constant term is just the explicit cosmological 
constant. Finally, as usual, p m follows the standard evo- 
lution equation 

p m + 3H(p m + p. m ) =0, (9) 

while p r (standard-model radiation) follows 

p r + 3H(p r +p r ) =0. (10) 

If we require expressions ([7]) to coincide with the stan- 
dard Friedmann equations, in units where c = 1 we set 
MM: 



In our work the running of A is not a problem, since 
the whole relevant cosmological history, that is after in- 
flation, is obviously inside the IR limit, that is with A = 1. 
On the other hand, the divergence of A from 1, that is 
the quantum gravitational features of the theory, become 
significant at very early times, that is close to the Big 
Bang or the cosmological bounce. This was additionally 
shown in (78|, where a detailed analysis based on SNIa, 
BAO and CMB observations, as well as BBN considera- 
tions, constrained A inside a narrow window around its 
infrared value: |A — 1| < 0.02. In summary, since in this 
work we are interested in post-inflation evolution, we set 
A = 1, and thus we simplify our notation using G for 
the (coincided gravitational and cosmological) Newton's 
constant. 



B. Beyond Detailed Balance 

The aforementioned formulation of Hofava-Lifshitz 
cosmology has been performed under the imposition of 
the detailed-balance condition. However, in the litera- 
ture there is a discussion whether this condition leads to 
reliable results or if it is able to reveal the full informa- 
tion of Hofava-Lifshitz gravity [U HH. Therefore, one 
needs to investigate also the Friedman equations in the 
case where detailed balance is relaxed. In such a case one 
can in general write @, ©, [TJ SI 51 : 



H 2 = 



2ct 



(3A - 1) 
2 

(3A - 1) 

(To 



11 

6 
A" 



Pr , 

°zK 2 
6a 4 



a 4 K 



6a 6 



3(3A- 1) a 2 



(13) 



16tt(3A - 1) 



_«VA_ 

8(3A- 1) 



(11) 



where G COS mo is the "cosmological" Newton's constant, 
that is the one that is read from the Friedmann equa- 
tions. We mention that in theories with Lorentz invari- 
ance breaking G COS mo does not coincide with the "grav- 
itational" Newton's constant G gra v, that is the one that 
is read from the action, unless Lorentz invariance is re- 
stored |l!3| . For completeness we mention that in our 
case 



Ggrav - 32ir , 



(12) 



as it can be straightforwardly read from the action ^ 
(ou r definitions of G COS mo, G grav coincide with those of 
[16l [Hj]). Thus, it becomes obvious that in the IR 
(A = 1), where Lorentz invariance is restored, G cosmo 
and G erav coincide. 



H 



-H 2 = - 



3(T 



(3A-1) 
3 



(3A-1) 

(T-2 



WraPrn + W r p r 
2 



6 



&3K <j 4 K 



18a 4 



6a 6 



A" 



6(3A- 1) a 2 



(14) 



where (To = k 2 /12, and the constants <Ji are arbitrary 
(with 02 being negative and 04 positive). As we observe, 
the effect of the detailed-balance relaxation is the decou- 
pling of the coefficients, together with the appearance of 
a term proportional to a~ 6 . 

Finally, if we force (|13 |) . (|T4l) to coincide with the stan- 
dard Friedmann equations, we obtain: 



6er 



G cosmo - ^ (3A _ x) 
a 2 = -3(3A - 1), 



(15) 



while in this case the "gravitational" Newton's constant 
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Gg ra v reads @: 



r - 6cr ° 
Ugrav " 16tt" 



(16) 



Similarly to the detailed balance case, in the IR (A = 
1) G cosmo and G grav coincide, and thus in the following 
we will use G to denote the (coincided gravitational and 
cosmological) Newton's constant. 



III. GENERALIZED CHAPLYGIN GAS IN 
HORAVA-LIFSHITZ COSMOLOGY 



Thus, since in this work we are interested in investigating 
the generalized Chaplygin gas in Hofava-Lifshitz cosmol- 
ogy, we first separate the total matter energy density p m , 
that was present in the Friedmann equations of the pre- 
vious section, into the baryonic pb and the dark matter 
Pdm parts, and then we use the generalized Chaplygin gas 
energy density p c instead of pdm- Therefore, using also 
the definitions (fTTj) . the Friedmann equations of 
the detailed balance case take the form: 



H 



inG 



\Pb 



Pr\ 



2Aa 4 



K 
~7p- 



(22) 



As discussed in the introduction, Chaplygin gas scenar- 
ios have been studied extensively as possible candidates 
for a unified description of dark matter and dark energy 
through a single fluid. The generalized Chaplygin gas is 
a fluid with an equation of state given by: 



A 

Pc = ? , 

Pc 



(17) 



with A a positive constant and < (3 < 1. p c and p c 
denote the energy density and pressure of the Chaplygin 
gas, and henceforth a subscript "c" will denote quantities 
pertaining to the Chaplygin gas. The simple Chaplygin 
gas is the special case of (3 = 1. Finally, considering the 
covariant conservation of energy-momentum, p c satisfies 
the standard evolution equation 



p c + 3H(p c + p c ) = 0. 



(18) 



Equations (JT7J) and ([18]) can be easily solved to yield 
MM: 



PcO 



A. 



1 - A„ 



a 3(l+/3) 



(19) 



where p c o denotes the present-day density of the Chaply- 
gin gas (in what follows, a subscript will always denote 
the present-day value of a quantity), and we have intro- 
duced the parameter 

A S ^A P ^\ (20) 
Thus, (fTTf provides the pressure of the Chaplygin gas as 

1-As 



Pc = - Pea A s 



A, 



a 3(l+/3) 



3 

' i+f) 



(21) 



Clearly, A s is therefore just the present value of the 
cquation-of-state parameter of the generalized Chaply- 
gin gas. Furthermore, j3 is related to the sound velocity 
of the C hapl ygin gas at the present time, which is given 
by (3A S [ToJ. 

The most important advantage of the Chaplygin gas 
scenario is that it behaves like dark matter at early times 
(a -C 1) and as a cosmological constant at late times 
(n» 1), smoothly interpolating between the two phases. 



H 



-4ttG 



1 




r k 2 


3A" 


Pc + gPr 




4Aa 4 


~ T 



K 
(23) 

Similarly, for the beyond-detailed-balance case the Fried- 
mann equations (fT3")) , (|14[) . using also the identifications 
(|T5l) . become 



H 2 = 



IttG 



[pb + Pc + Pr. 

2 IVi 
~6~ 



(3A-1) 



<j 3 K 2 ^ <r 4 K 
6a 4 a 



H 



-H 



-AttG 



(3A - 1) 



11 
6 



1 

a 3 K 2 
18a 4 



6a 6 



K 
7a 



(24) 



6a 6 



K 
2a 2 



.(25) 



IV. OBSERVATIONAL CONSTRAINTS 

In the conventional gravitational background of gen- 
eral relativity, the cosmological scenarios of generalized 
Chaplygin gas have been extensivel y constrained using 
obse rvational data. In particular, in jlOS . 1041 1071 . 1114 
Ill5| SNIa observations wer e used, in 107l Il08| data 
fro m la rge scale structure, in [lll[Il7| CMB data, while 
in [ll4| the authors used gravitational lensing observa- 
tions. Finally, constra ints from com bined data sources 
have been obtained in |l!8| and |l!9| . 

In the present work we will consider observations in or- 
der to constraint the scenario of a generalized Chaplygin 
gas in a universe governed by Hofava-Lifshitz gravity. 
In particular we will use data from SNIa, BAO, CMB 
and BBN to constrain the Chaplygin gas parameters to- 
gether with the parameters of Hofava-Lifshitz gravity. 
For completeness, we perform our analysis separately 
for the detailed-balance and the beyond-detailed-balance 
scenario. 

A. Constraints on the Detailed-Balance scenario 

We first consider the detailed-balance version of the 
theory. It proves more convenient to use the redshift z 
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instead of the scale factor as the independent variable, 
through the relation 1 + z = do/a = 1/a. Furthermore, 
we use the dimensionless density parameters 



o - 8nG o K 

si) = n Pi, liK = —- 



3H2 rt> - F2a2 , 

and we introduce the dimensionless parameter 

A 



2HI 



Finally, we use the dimensionless expansion rate 

_H(z) 



E{z) 



(26) 



(27) 



(28) 



Using the above definitions, the Friedman equation 
221) becomes: 



e 2 {z) = n M (i + z) 3 + n c0 F( z ) + n r0 (i + z) A 



+ n K0 (i + z) 2 + 



o 2 

^(l + z) 4 



where 



F(z)= A s + (l-A s )(l + z) 



3(1+0) 



(29) 



(30) 



Obviously, (|2"9")l is subject to the present-time constraint 
i?(z = 0) = 1, which leads to the condition: 



slfcO + il c Q + il r o + + W + 



O 2 
4w 



(31) 



As we have already mentioned above, the term 
f2 2 ^ /(4w) is the coefficient of the dark radiation term, 
which is a characteristic feature of the Hofava-Lifshitz 
gravitational background. Since this dark radiation com- 
ponent has been present also during the time of nucle- 
osynthesis, it is subject to bounds from Big B ang Nucle- 
osynthesis (BBN). As discussed in detail in [73], if the 
upper limit on the total amount of Hofava-Lifshitz dark 
radiation and kination-l ike (a qu intessence field domi- 
nated by kinetic energy 12dl Il2ll | ) components allowed 
during BBN is expressed thro ugh the p arameter AN V of 
the effective neutrino species 
the following constraint: 



122T4T25 1 , then we obtain 



O 2 
4w 



0.135AA^Sl 



rO- 



In this work, in order to ensure consistency wit h BB N, 
we adopt an upper limit of AN V < 2.0 following jl24l |. 

In most studies o f dark ene rgy models it is usual to ig- 
nore curvature (e.g. |l29l4l37| ). especially concerning ob- 
servational constraints. This consideration is well mo- 
tivated since most inflationary scenarios predict a high 
degree of spatial flatness, and moreover the CMB data 
impose stringent constraints on spatial flatness in the 



context of constant-w models (fo r ex ample a combina- 
tion of WMAP+BAO+SNIa data PH provides the tight 
simultaneous constraints —0.0179 < FIko < 0.0081 and 
-0.12 < 1 + w < 0.14, both at 95% confidence). 

However, it is important to keep in mind that due 
to degeneracies in the CMB power spectrum (see |l39j 
and references therein), the limits on curvature depend 
on assumptions regarding the underlying dark energy 
scenario. For example, if we use a linear w (that is 
w (a) = wq + (1 — a) w a ) instead of a constant one, the 
error on f^ p is of th e ord er of a f ew percent, that is 
much larger [140Hl4^ (see [IZMI^ for the constraints 
on curv ature for different parameterizations). Addition- 
ally, in 145j it was shown that for some models of dark 
energy the constraint on the curvature is at the level of 
5% around a flat universe, whereas for others the data 
are c onsistent with an open universe with Qko ~ 0.2. In 



14l| it was proposed that geometrical tests, such as the 



combination of the Hubble parameter H(z) and the an- 
gular diameter distance Da{z), using (future) data up to 
sufficiently high redshifts z ~ 2, might be able to disen- 
tangle curvature from dark energy evolution, t hough not 
in a model-independent way. Furthermore, in [l46l Il47| 
the authors highlighted the pitfalls arising from ignor- 
ing curvature in studies of dynamical dark energy, and 
recommended to treat Qko as a free parameter to be 
fitted along with the other model parameters. Lastly, 
note that in the present work the spatial curvature plays 
a very crucial role, since Hofava-Lifshitz cosmology co- 
incides completely with ACDM if one ignores curvature 
[il], HI]- Therefore, and following the discussion above, 
we prefer to treat £Iko as a free parameter. 

In summary, the scenario at hand involves the follow- 
ing parameters: flbo, fl r o, sl c o, sl^o, w, AN„, Hq, A s 
and /3. We fix the parameters sl m o(= ^bo + sl c o), f^n , 
H and O r0 at their 7-year WMAP best-fit values [L48j . 
namely sl m0 = 0.27, sl r0 = 8.14 x 10~ 5 and H =71.4 
Km/sec/Mpc. Therefore, there are five remaining pa- 
rameters, namely flko, w > AN V , A s and j3, which are 
subject to the two constraint equations (f3Tj) and (|32|) . 
We choose to treat AN Vi A s and /3 as free parameters 
and use these constraints to eliminate &ko and w for a 
given choice of curvature, as: 

w (K; AN„,A s ,f3) = 1 - Q m0 - (1 - .135A^)sl r0 

-0.73 sgn(A-) y/AN^y/ fi r0 - - sl 2 (33) 

and 



(32) \n K0 (AN v ,A s ,l3)\ = ^0.b4AN„n rO u(K;AN„,A s ,p). 

(34) 

Note that in order to obtain w and £Iko from the above 
equations (|3"3"|) and (f3~4"|). the type of curvature K has to 
be chosen a priori. Therefore, we treat positive {K > 0) 
and negative (K < 0) curvature separately. Finally, note 
that according to (|34[) if AJV„ = then it forces the 
spatial curvature to be zero. 

Let us now proceed to constrain the three free param- 
eters AN„, A s and j3. We perform a likelihood analysis 



6 





0.01 0.02 0.03 0.04 0.05 0.06 0.07 
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FIG. 1: (Color online) Detailed Balance: Likelihood contours of A s and AN V for different choices of f3 as discussed in the 
text. The lightest region is excluded at the 2o level, and the darker region is excluded at the la level. The darkest region is not 
excluded at either confidence level. 



using the SNIa, BAO and CMB data (see the appendix 
of [73] for the details of the procedure) . In Fig. [T] we 
present the likelihood contours of A s and A7V„ for three 
different choices of j3, one at its best-fit value, and the 
other two at the extremes of its allowed values. As we ob- 
serve, the scenario at hand is in general compatible with 
observations. Furthermore, note that the contours are 
smaller for negative curvature, indicating that negative 
curvature is generally disfavored for these models. Fi- 
nally, the best-fit values together with the corresponding 
la bounds arising from the analysis, for the two choices 
of curvature, are shown in Table HI 





K < 


K > 


A s 


1.96 x 1(T 13 ; (0,0.04) 


3.75 x 10- 11 ; (0,0.04) 


AN,, 


0; (0,0.09) 


2.75 x 10~ 2 ; (0,0.09) 


P 


3.21 x 10" 5 ; (0,1) 


3.89 x 10" 3 ; (0,1) 



TABLE I: Best-fit values and \o confidence intervals for the 
free parameters A s , AN„ and f3 in the detailed balance sce- 
nario, arising from a likelihood analysis using SNIa, BAO and 
CMB data. 



B. Constraints on the Beyond-Detailed-Balance 
scenario 



In this subsection we consider the cosmology of a gen- 
eralized Chaplygin gas in Hofava-Lifshitz gravity where 
the detailed balance condition has been relaxed. Our cos- 
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mic fluid consists of the baryonic matter, the standard- 
model radiation and the Chaplygin gas. Starting with 
the Friedmann equation (fT3|) and proceeding similarly to 
the previous subsection, we can write the dimensionless 
Hubble parameter E(z) as: 



e 2 (z) =n w (i + zf + n cQ F( z ) 
+ n r0 (i + zf + n K0 (i + zf 

+ wi + w 3 (l + zf + w 4 (l + zf 
where as usual 

F(z)= \A s + (l-A s )(l + zf^ 



(35) 



(36) 



In ()35[) we have introduced the dimensionless parameters 
Wi, and W4, related to the model parameters ax, 03 
and (T4 through: 



CJ4 



6^o 2 



6 

6 ' 



(37) 



Furthermore, we consider the combination lo^ to be posi- 
tive, in order to ensure that the Hubble parameter is real 
for all redshifts (note that 0J4 > is required also for the 
stability of the gravitational perturbations of the theory 
0EI)- For convenience we moreover assume cr 3 > 0, 
that is uj 3 > 0. 



The present scenario involves the following parameters: 
H , O 60 , ^co, ^ro, ^ko, "1, co 3 and w 4 , A s and /3. Simi- 
larly to the detailed-balance case, we fix the parameters 
il m n ,flhq , Ho and Sl r o a t their 7-year WMAP best-fit val- 
ues [148l |. Thus, the remaining free parameters are Qko, 
ojx, W3 and LU4, A s and /3, which are subject to the same 
two constraints as discussed in the previous subsection. 
The first one arises from the Friedman equation at z = 0, 
which leads to 



Q m o + fir-o + Qko + wi + W3 + W4 = 1. 



(38) 



We use this constraint to eliminate the parameter wx- 

The second constraint arises from BBN considera- 
tions. The term involving W3 represents the usual dark- 
radiation component. In addition, the W4-term represents 
a kination-like compone nt (a qui ntessence field domi- 
nated by kinetic energy |l2dl . Il2l| ). also called a "stiff 
fluid". If AN V represents the BBN upper limit on the 
total energy density of the universe beyond standard 
model constituents, then as shown in the Appendix of 



[77j . we hav e the follo wing constraint at the time of BBN 
(z = zbbn) [Mill: 



uj 3 + u± (1 + z B bn) = w 3max = 0.135AA^ 



r0- 



(39) 



It is clear that BBN imposes an extremely strong con- 
straint on CJ4, since its largest possible value (correspond- 
ing to W3 = 0) is ~ 10 -24 . This is in agreement with 
precision constra ints on stiff fluid densities at the time of 
BBN derived in [l26l |. Finally, W3 max denotes the upper 
limit on W3. In the following, we use expression (|39p to 
eliminate W4. For convenience, instead of w 3 we define a 
new parameter 



cj 3 



W3n 



(40) 



which has the interesting physical meaning of denoting 
the ratio of the energy density of the Hof ava-Lifshitz dark 
radiation to the total energy density of Hofava-Lifshitz 
dark radiation and kination-like components at the time 
of BBN. 

In summary, in the scenario at hand we have the fol- 
lowing five free parameters: AN U , Qko, a, A s and /3. We 
use SNIa, BAO and CMB data to perform a likelihood 
analysis, and we construct the likelihood contours for the 
parameters A s and AiV„ for various fixed choices of the 
other parameters a, /3 and £Iko- 

In Fig. [2] we present the contours of A s and AN V 
obtained by varying a around its best-fit value, while 
keeping /3 and Qko fixed at their best-fit values. We 
find that as a increases, the contours shrink and rotate 
counter-clockwise, placing tighter constraints on AN V . 

In Fig. [3] we present the contours of A s and AN U for 
four different choices of spatial curvature, with a and (3 
fixed at their best-fit values. We observe that curvature 
constraints arc quite tight. For example the 1-a regions 
shrink as |fifeo| increases from 0.001 to 0.002 (as shown 
in the figure) and they almost completely disappear for 
|fifco| ^ 0.003 (which is not shown). 

In Fig. [4] we show the likelihood contours of A s and 
A7V„, with a and £Iko fixed to their best-fit values, and 
two extreme choices of (3. We deduce that the contours 
expand slightly for the larger value of (3. 

Finally, the best-fit values, along with the correspond- 
ing la confidence intervals are presented in Table [XXJ 



A a 


2.47 x 1(T 10 ; (0,2) 


AN, 


5.23 x 1(T 4 ; (0,2.0) 


a 


5.17 x 10" 4 ; (0,1) 


Q,K0 


1.29 x 10~ 3 ; (-0.003,0.003) 


P 


4.89 x 10" 5 ; (0,1) 



TABLE II: Best-fit values and la confidence intervals for the 
free parameters A s , AN V , a, Qko and /3, in the beyond- 
detailed-balance scenario, arising from a likelihood analysis 
using SNIa, BAO and CMB data. 
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FIG. 3: (Color online) Beyond Detailed Balance: Likelihood contours of A s and AN U for four choices ofQfco, with the other 
parameters fixed to their best-fit values. Color scheme as in Fig. QJ 
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total pressure and energy density given by 



FIG. 4: (Color online) Beyond Detailed Balance: Like- 
lihood contours of A s and AN V for two extreme different 
choices of ft, with the other parameters fixed to their best-fit 
values. Color scheme as in Fig. QJ 



V. COSMOLOGICAL IMPLICATIONS 



In the previous sections we constructed cosmological 
scenarios of a generalized Chaplygin gas in a universe 
governed by Hofava-Lifshitz, and we used observational 
data in order to impose constraints on the model pa- 
rameters. In recent past, people have discussed different 
cosm ologi cal scenarios in Hofava-Lifshitz gravity. Maeda 
et al jl27l | have recently studied the bouncing as well as 
oscillating universe in Hofava-Lifshitz gravity in presence 
of both curvature and cosmological constant. One very 
interesting feature of their study is the possibility of a 
quantum tunneling from the oscillating spacetime to an 
inflationary sce nario . In another interesting work, Berto- 
lami and Zarro jl28l | have discussed quantum cosmology 
via minisupcrspacc model for projectablc Hofava-Lifshitz 
gravity without detailed balanced condition. They con- 
sidered the presence of a cosmological constant and when 
it is positive they recovered the classical GR solution with 
accelerating expansion. In the present section we discuss 
some of the cosmological implications in our case where 
we consider the presence of generalised Chaplygin gas . 
In particular, we focus on the evolution of the dimension- 
less expansion rate defined in (|2"8|) . and on the evolution 
of the equation-of-state parameter of the total cosmic 
fluid of the universe, defined as w = Ptot/ptot, with the 



Plot 



1 2 



Plot = Pc+ Pb+ Pr + 



K 2 
Aa 4 " 

m 2 



3A 



Aa 4 



3A 



for the detailed-balanced scenario, and 
1 



Plot — Pc + 



3 pr 



Plot = Pc 



Pb 



6ct 



a 3 K 2 


<7 4 K 


18a a 4 


6<7o a6 


a 3 K 2 




OfToa 4 


6<7o a6 



(41) 



(42) 



(43) 



(44) 



for the beyond-detailed-balance one, as it can be eas- 
ily extracted from the corresponding Friedmann equa- 
tions. Therefore, replacing the scale factor by the red- 
shift, and using the density parameters (|26[) at present, 
as well as the Hubble parameter as a function of the red- 
shift (from relations (|29p and (|35[) . for the detailed and 
beyond detailed-balance case respectively), we straight- 
forwardly acquire the total equation-of-state parameter 
function of the redshift w(z). 



A. Detailed Balance 

Let us first examine the evolution of the total equation- 
of-state parameter w{z) of the universe. In order to make 
our presentation simpler we will use the best-fit value set 
of Table fl] as our fiducial choice for the various model 
parameters, denoting by Wbf(z) the corresponding w{z). 
That is, the subscript "bf" for a parameter will imply 
the best-fit value of the parameter and for a variable will 
imply that all parameters determining that variable have 
been fixed to their best-fit values. In Fig. [5] we present 
the evolution of Wbf(z). 



As we observe, at very early times the cosmic equation 
of state is close to 1/3, since the standard-model radi- 
ation and the dark radiation from the Hofava-Lifshitz 
gravitational background dominate, as expected. At in- 
termediate redshifts, the Chaplygin gas dominates and it 
behaves like a dust for quite a long time, leading to an 
equation of state close to zero. Finally at late times, the 
de Sitter phase of the Chaplygin gas dominates. Inter- 
estingly, this evolution is consistent with the thermal his- 
tory of our universe, with the succession of the radiation, 
matter, and "dark energy" epochs, and this feature acts 
as an advantage of the scenario at hand. Furthermore, 
note that we do not need any additional mechanism to 
describe the late-time universe acceleration (that is the 
dark energy) , since this is obtained through a unified way 
by the generalized Chaplygin gas. That is, this crucial 
property of the Chaplygin gas is still exhibited by the 
present Hofava-Lifshitz embedded scenario. 
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FIG. 5: (Color online) Detailed Balance: Evolution of the 
fiducial total equation- of- state parameter of the cosmic fluid 
Wbf(z), for positive and negative curvatures, with all parame- 
ters fixed to best-fit values presented in Tabled 





FIG. 6: (Color online) Detailed Balance: Evolution of the 
fiducial dimensionless Hubble parameter Ebf(z) for positive 
and negative curvatures with all parameters fixed to best-fit 
values presented in Tabled 



We now focus on the evolution of the dimensionless 
expansion rate E(z) denned in In order to acquire 

the basic information, in Fig. [5] we depict the evolution 
of the fiducial expansion rate Ehf(z), that is with all pa- 
rameters fixed to their best-fit values presented in Table 
HJ As expected, Ebf(z) starts from very large values, and 
it steadily decreases towards its present value 1. 

It would be interesting to investigate how the expan- 
sion rate is affected by each parameter. In order to quan- 
tify this, we define the quantity Q(z), which is the frac- 
tional change in the expansion rate Ehf(z) as each pa- 
rameter is varied within its allowed limits, while all the 
other parameters are kept fixed at their best-fit values, 
that is 



E(z)-E bf (z) 
Q\ z > = W~f7\ ■ 

Ebf{z) 



(45) 



Note that A s and f3 have well-defined upper and lower 
limits, while for AN V we use the standard BBN upper 
limit of AN V < 2.0 [l22|. In Fig. we present Q(z) 
for all the parameters of the scenario. As we see, both 
A s and AN U can have a strong impact on the expansion 
rate, while j3 has a very weak impact. The fact that the 
effect of (3 is small was expected, since in the contour 
plots of Fig. [T] we did not obtain a significant change 
upon varying f3. 

Finally, from Fig. [7] we deduce that the impact of 
both the generalized Chaplygin gas parameters A s and 
/3 is realized mainly at low redshifts, while it weakens at 
high ones. This is explained by the fact that at high red- 
shifts the Chaplygin gas is subdominant and thus the ex- 
pansion dynamics is driven by the (standard-model and 
dark) radiation. On the other hand, AN V , which is di- 
rectly related to the amount of dark radiation, can have 
a significant impact on the expansion rate at arbitrarily 
high redshifts. 



B. Beyond Detailed Balance 

Let us now repeat the analysis of the previous subsec- 
tion, in the beyond-detailed-balance version of the exam- 
ined scenario. 

In Fig. |8]we depict the evolution of Whf(z), that is the 
total equation-of-state parameter of the universe, with 
all the model parameters fixed to their best-fit values 
presented in Table HU As we observe, the behavior of 
Wbf(z) is similar to that of the detailed-balance case, that 
is we obtain an initial radiation era, a relatively long 
matter one, and finally an era of accelerated expansion. 

Concerning the dimensionless expansion rate of the 
universe E(z), in Fig. |9]we depict its fiducial value Ebf(z) 
as a function of the rcdshift. Again we see that Ebf(z) 
starts from very large values and it steadily decreases, 
compatibly with observations. 
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FIG. 8: (Color online) Beyond Detailed Balance: Evo- 
lution of the fiducial total equation- of -state parameter of the 
cosmic fluid Wbf(z), with all parameters fixed to best-fit values 
presented in Table [771 



FIG. 9: (Color online) Beyond Detailed Balance: Evolu- 
tion of the fiducial dimensionless Hubble parameter Ebf(z), 
with all parameters fixed to their best-fit values presented in 
Tabled 
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FIG. 10: (Color online) Beyond Detailed Balance: Frac- 
tional variation Q(z) of the dimensionless Hubble parameter 
Ebf(z), upon varying the parameter A s within its allowed lim- 
its, as indicated in the graph. All the other parameters are 
kept fixed to their fiducial values presented in TableWl 
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FIG. 11: (Color online) Beyond Detailed Balance: Frac- 
tional variation Q(z) of the dimensionless Hubble parameter 
Ebf(z), upon varying the parameter AiV„ within its allowed 
limits, as indicated in the graph. All the other parameters are 
kept fixed to their fiducial values presented in TableUK 



We now proceed to the investigation of the effect that 
the various model parameters have on the expansion rate. 
In Figures HOITmI wc demonstrate the fractional change 
Q{z) of the expansion rate E b f(z) as each parameter is 
varied within its allowed limits, as it is indicated in the 
graphs, while all other parameters are kept fixed at their 
best-fit values. From Fig. [Tniwe deduce that A s can have 
a strong impact on the expansion rate at low redshifts, 
while from Fig. [TT] we see that AiV„ can affect signifi- 
cantly the expansion rate at large redshifts. On the other 
hand, form Figures [T^J [T3] and we observe that f^o 
and especially the parameters a and f) have negligible 
impacts. Similarly to the detailed balance case, the fact 
that the effect of f3 is small was expected, since in the 
contour plots of Fig. |4] we did not obtain a significant 
change upon varying (3. 
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FIG. 13: (Color online) Beyond Detailed Balance: 

Fractional variation in the dimensionless Hubble parameter 
Ebf(z), upon varying the parameter a within its allowed lim- 
its, as indicated in the graph. All the other parameters are 
kept fixed to their fiducial values presented in TableUH 
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FIG. 12: (Color online) Beyond Detailed Balance: Frac- 
tional variation Q(z) of the dimensionless Hubble parameter 
Ebf(z), upon varying the parameter Q,ko as indicated in the 
graph. All the other parameters are kept fixed to their fiducial 
values presented in TableUH 
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FIG. 14: (Color online) Beyond Detailed Balance: 

Fractional variation in the dimensionless Hubble parameter 
Etf(z), upon varying the parameter ft within its allowed lim- 
its indicated, as in the graph. All the other parameters are 
kept fixed to their fiducial values presented in Table\U\ 



The explanation of the effect of the various parameters 
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on the expansion rate is easily obtained, knowing the 
behavior of the various terms as a function of the redshift. 
In particular, at high rcdshifts the radiation (standard- 
model and dark one) dominates and drives the expansion, 
while the Chaplygin gas is subdominant. However AN U , 
which quantifies the Hofava-Lifshitz dark radiation and 
kination-like components, can have a significant impact 
on the expansion rate at arbitrarily high redshifts. 



VI. CONCLUSIONS 

In this work we investigated the cosmological scenario 
of generalized Chaplygin Gas in a universe governed by 
Hofava-Lifshitz gravity, both in the detailed-balance as 
well as in the beyond-detailed balance version of the the- 
ory. Furthermore, in order to obtain a realistic cosmol- 
ogy we have included the baryonic matter and standard- 
model radiation sectors. 

After extracting the cosmological equations, we used 
data from SNIa, BAO and CMB observations, as well 
as arguments from Big Bang Nucleosynthesis, in order 
to impose constraints on the various parameters of the 
scenario. For the detailed-balance case we have three 
free parameters, namely the exponent f3 of the gener- 
alized Chaplygin Gas, its present equation-of-state pa- 
rameter value A s , and the effective neutrino parame- 
ter AA„ which quantifies the total amount of Hofava- 
Lifshitz dark-radiation and kination-like components al- 
lowed during BBN. The best fit values of the parame- 
ters together with the corresponding lu confidence in- 
tervals, arisen from the likelihood analysis for open and 
closed universe, were shown in Table HI while in Fig. Q] 
we presented the corresponding likelihood contours. We 
deduced that the scenario at hand is compatible with 
observations, however the data lead to strong bounds on 
A s and AN V and the spatial curvature £Iko- Finally, the 
scenario at hand depends only slightly on /3. 

In the beyond-detailed-balance scenario the free pa- 
rameters are five, namely (3, A s , AN V , the curvature 
density Qjfo, an d the parameter a, which is the ratio 
of the dark-radiation energy density to the sum of dark- 
radiation and kination-like energy densities at the time of 
BBN. The best fit values of the parameters together with 
the corresponding la confidence intervals, arisen from 
the likelihood analysis, were shown in Table [TTJ Addi- 
tionally, in Figures [2H] we presented the likelihood con- 
tours of the model parameters. From this analysis we 
deduced that the entire (consistently with BBN) range 
of < AA„ < 2.0 is allowed, for suitable choices of the 
other parameters. We also found tight constraints on the 
curvature and on A s . 

After the observational elaboration, we investigated 
the cosmological implications of the examined scenar- 
ios. In particular, we focused on the evolution of the 
total equation-of-state parameter of the universe, and of 
the expansion rate. In the detailed-balance case, we saw 
that in general at very early-times is radiation that domi- 



nates, at intermediate redshifts is the Chaplygin gas that 
dominates, behaving like matter for a long time, and fi- 
nally at late times the dominant Chaplygin gas behaves 
like dark energy, triggering the accelerating expansion 
(see Fig [SJ. This evolution is consistent with the thermal 
history of the universe, and this is an advantage of the 
present scenario. Note moreover that the present accel- 
erated era is described in a unified way, without the need 
of any additional mechanism. Concerning the qualitative 
dependence on the various parameters, we found that the 
expansion rate has a strong dependence on both A s and 
AA„ at low redshifts, while at high redshifts, where the 
Chaplygin gas is not dominant, this dependence weakens 
and disappears (see Fig [7]). Lastly, the dependence on (3 
is very weak at all redshifts. 

In the beyond-detailed-balance scenario we found a 
similar behavior for the total equation-of-state param- 
eter of the universe, that is we obtained successively a 
radiation, a matter, and a dark energy era (see Fig [5]). 
Qualitatively, A s has a weak impact on the expansion 
rate at lower redshifts, AA„ has a significant effect at 
large redshifts, while £Iko and especially a and (3 have 
a negligible impact on the expansion rate at all redshifts 
(see Figures [T(IT4|) . 

In summary, the generalized Chaplygin gas in Hofava- 
Lifshitz gravitational background is compatible with ob- 
servations, and can successfully reproduce the expansion 
history of the universe. However, we should mention 
that the present analysis does not enlighten the discus- 
sion about the possible conceptual problems and insta- 
bilities of Hofava-Lifshitz gravity, nor it can address the 
questions concerning the validity of its theoretical back- 
ground, which is the subject of interest of other studies. 
It just analyzes the phenomenological consequences and 
the cosmological implications of the generalized Chaply- 
gin Gas in such a gravitational background, and thus its 
results can be taken into account only if Hofava-Lifshitz 
gravity passes successfully the necessary theoretical tests. 

In the same lines, the present work does not address 
the problem of undesirable instabilities and oscillations 
in the matter power spectra, that GCG scenario faces in 
standard Einstein gravity. Although, as discussed in the 
Introduction, there are several approaches which try to 
solve this problem in conventional gravity, it may be quite 
interesting to see whether it can be addressed in the con- 
text of Hofava-Lifshitz gravity. Let us make som e com- 
ments in this perspective. As discussed in |98lll03l |. GCG 
can be modelled in terms of scalar fields having both 
canonical as well as non-canonical kinetic energy terms, 
and thus perturbing the GCG one has to essentially per- 
turb these scalar fields. In the case of Hofava-Lifshitz 
gravity one expects in general the linear perturbations 
to be quite different due to higher-order curvature-terms 
present in the action. These higher-order gradient terms 
can in principle generate non-adiabatic pressure pertur- 
bation, which can help to cure the instabilities or os- 
cillations in the power spectra even in the large scales. 
The effect can be prominent with scalar fields having 
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non-canonical kinetic energy terms, for instance of Dirac- 
Born-Infcld form, which is typical for GCG-likc equation 
of state. 

In conclusion, the generalized Chaplygin gas scenario 
in Hofava-Lifshitz cosmology, can have rich cosmologi- 
cal consequences. The present study is a first step in 
this direction, where we have considered the background 
evolution and have confronted the model with a variety 
of currently available observational data. The next step 
is to consider the inhomogeneous GCG model in Hofava- 
Lifshitz gravity, which will be our goal in the near future. 
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